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Abstract 

Background: Hepatocellular carcinoma is a primary tumor of the liver and involves different treatment modalities 
according to the tumor stage. After local therapies, the tumor evaluation is based on the mRECIST criteria, which 
involves the measurement of the maximum diameter of the viable lesion. This paper describes a computed 
methodology to measure through the contrasted area of the lesions the maximum diameter of the tumor by a 
computational algorithm. 

Methods: 63 computed tomography (CT) slices from 23 patients were assessed. Non-contrasted liver and HCC typical 
nodules were evaluated, and a virtual phantom was developed for this purpose. Optimization of the algorithm detection 
and quantification was made using the virtual phantom. After that, we compared the algorithm findings of maximum 
diameter of the target lesions against radiologist measures. 

Results: Computed results of the maximum diameter are in good agreement with the results obtained by radiologist 
evaluation, indicating that the algorithm was able to detect properly the tumor limits. A comparison of the estimated 
maximum diameter by radiologist versus the algorithm revealed differences on the order of 0.25 cm for large-sized 
tumors (diameter > 5 cm), whereas agreement lesser than 1.0 cm was found for small-sized tumors. 

Conclusions: Differences between algorithm and radiologist measures were accurate for small-sized tumors with a 
trend to a small decrease for tumors greater than 5 cm. Therefore, traditional methods for measuring lesion diameter 
should be complemented non-subjective measurement methods, which would allow a more correct evaluation of 
the contrast-enhanced areas of HCC according to the mRECIST criteria. 



Background 

Hepatocellular carcinoma (HCC) is the fifth most com- 
mon cancer in men and the seventh most common in 
women [1]. Its incidence is highest in regions where 
hepatitis B virus is endemic [1]. In the United States, 
deaths caused by hepatitis C virus (HCV) -related HCC 
are rapidly rising. In the past two decades, the incidence 
of HCC in the United States has tripled, but the 5-year 
survival rate has remained below 12% [2]. The greatest 
proportional increase in cases of HCC has been seen 



* Correspondence: matheus@ibb.unesp.br 

'institute de Biociencias de Botucatu, Departamento de Fisica e Biofisica, UNESP - 
Universidade Estadual Paulista, Distrito de Rubiao Junior S/N, Botucatu, 18618-000 
Sao Paulo, Brazil 

Full list of author information is available at the end of the article 



among Hispanics and whites between 45 and 60 years 
old [2]. 

In general, HCC diagnosis is based on noninvasive im- 
aging tests [3,4]. In patients with cirrhosis and a focal 
hepatic lesion > 2 cm, the diagnosis may be confidently 
established on the basis of typical imaging features showing 
areas of arterial enhancement and regions promptly 
"washed out" (fainter than the liver tissue) in the venous 
or delayed phase of four-phase multidetector computed 
tomography (CT) exam (where the four phases are unen- 
hanced, arterial, venous, and delayed) [3,5]. 

Orthotopic liver transplantation (OLT) is the recom- 
mended treatment when the tumor is within specific 
criteria, that based on the tumor maximum diameters 
and the absence of tumor spread outside the liver [3]. 
Previous studies support the effectiveness of OLT in 
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patients meeting the Milan criteria adopted by the United 
Network for Organ Sharing (UNOS) [4,5]. The Milan cri- 
teria states that orthotopic liver transplantation is recom- 
mended only if the patient have a solitary HCC nodule 
with a diameter < 5 cm or no more than 3 nodules with 
diameters < 3 cm [5,6]. 

When any kind of local therapy is used in order to re- 
duce the tumor size, the HCC nodules are evaluated 
according to the mRECIST criteria, which require that 
only well-delineated, arterially enhanced lesions could 
be selected as viable tumor tissue [7]. Using Computed 
Tomography (CT) images, the viable tumor is measured 
in the arterial phase, with highest distinction between 
the viable vascularized tumor nodule and necrotic tisue 
(non enhanced region) [8-11]. The longest diameter of 
the viable tumor should be carefully assessed on the CT 
examination, since it can change the decision of which 
treatment will be further indicated to each patient [7]. 
The measurement of the viable tumor diameter should 
not include any major intervening areas of necrosis [7]. 
However, large nodules and/or tumors previously sub- 
mitted to local therapies are often filled by necrotic 
areas, making difficult to calculate the viable area of 
these tumors [8]. 

Before a curative treatment, many patients with HCC 
receive non-curative therapies in order to achieve the 
best survival rates, and one of the most used of them is 
the transarterial chemoembolization (TACE). This kind 
of treatment consists of infusions of a chemotherapeutic 
agent through the vessel that is the responsible by the 
tumor nutrition [9]. As soon the infusion is finished, the 
vessel is closed in order to isolate the chemotherapeutic 
agent within the tumor [9]. As a result, it causes a local 
destruction of the initial lesion, which can be seen as a 
necrotic area within the tumor. 

The TACE efficacy can be measured by the decrease of 
the target lesion, and when the tumor has enough reduc- 
tion. The patients can be considered at low disease stage, 
allowing them to receive OLT as a curative treatment 
[6]. However, there are some difficult to distinguish be- 
tween patients which are good candidates to OLT and 
those still needing to receive other sessions of TACE or 
another therapy, taking only into account a single meas- 
ure of the maximum tumor diameter. Other criterion, 
as the WHO criterion is based on the sum of bidi- 
mensional perpendicular products. Another one is the 
RECIST criterion which is based on the unidimesional 
quantity, the sum of the longest found diameters. Both 
criteria were designed for the evaluation of cytotoxic 
agents and not for local therapies. TACE induces tumor 
necrosis with or without changes in tumor size [12]. Be- 
cause of this, the WHO and RECIST criteria have been 
considered as suboptimal methods for tumor response 
assessment in HCC patients undergoing TACE. As a 



result, the European Association for the Study of the 
Liver (EASL) and the American Association for the Study 
of Liver Disease (AASLD) have proposed new methods, 
including the concept of viable enhancing lesion modify- 
ing WHO (EASL) and RECIST (mRECIST) criteria, re- 
spectively [8]. The EASL and mRECIST criteria resulted 
in a higher objective response rate and provided more 
reliable prognostic information, including survival, than 
conventional WHO and RECIST criteria. So, the mea- 
surements of the viable areas of HCC nodules are now 
the best way to access the TACE efficiency [9] . 

In this work, we measured the maximum diameter of 
the viable HCC by a computed algorithm in patient after 
TACE. For this purpose, an algorithm based on the discrete 
wavelet transform (DWT) was developed to quantify the 
enhanced areas of the liver on CT images using a non sub- 
jective way. 

Methods 

CT data sample 

The CT data for this study comprises 23 enhanced high- 
resolution abdomen CT exams offered by the Clinical 
Hospital of Botucatu Medical School. Ethical approval 
was granted by Research Ethical Commission of this in- 
stitute under the protocol number 485/2012. Written in- 
formed consent was obtained from the patient for the 
publication of this report and any accompanying images. 
These exams included CT scans performed after transar- 
terial chemoembolization of 23 patients who had not 
been treated by another kind of therapy. The hospital 
radiologist selected 63 slices which contains at least one 
contrasted part of a typical HCC lesion. There were se- 
lected the slices that had the highest amount of tumoral 
contrast-enhanced tissue. Each CT image measured 
512 x 512 pixels, the pixel size ranged from 0.73 to 
0.89 mm, with mean pixel value 0.80 mm. The HCC max- 
imum contrasted diameter size evaluated by an expert 
radiologist, with more than three years of experience, 
varies from 1.4 cm to 13.3 cm, central value of 4.7 cm and 
mean ± standarddeviation of 5.3 ± 3.0 cm. 

Patients 

We included patients undergoing 4-phases tomography at 
Botucatu Medical School, who have all of the following 
conditions: (i) more than 18 years old; (ii) undoubted diag- 
nosis of cirrhotic liver and HCC superior to 1 cm of diam- 
eter; (iii) at least one contrast-enhanced lesion at the 
arterial phase; (iv) rapid washout of the lesion during the 
venous phase; (v) without previous lipiodol use in HCC le- 
sions treated through transarterial chemoembolization. 

Highlighting structures using wavelets 

Multiscale contrast enhancement was used to highlight 
the tumor inside the liver. Similar to Fourier analysis, 
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the WT corresponds to a decomposition of the functional 
representation of the digital image. Whereas Fourier trans- 
formation represents the signal in the frequency domain, 
the WT provides a spatio-frequency decomposition of the 
signal [10,11]. 

This section describes image decomposition by multi- 
resolution analysis (MRA), which has the ability to sep- 
arate the decomposition into higher frequency bands 
and residuals (low-frequency components). This method 
offers multiresolution properties and highlights the char- 
acteristics of interest in the image. By introducing a high 
band-pass function ij/ and a low-pass scaling function <p, 
a one-dimensional (ID) signal J{x) can be decomposed 
by MRA as: 

/ 

/(*) = E E + E c /(*) C 1 ) 

k /=1 k 

where the first and second terms on the right-hand side 
of equation (1) represent the decomposed high- and 
low-frequency components of the image, respectively 
[13]. The wavelet coefficients eL(k) are given by the sca- 
lar products of the original image with the xf/j^ basis ele- 
ments in the pixel position k (covering the row image) 
and the chosen decomposition level / Q(/c) are the scal- 
ing coefficients, and <pj,k{x) are the respective scaling 
functions [12]. 

In the case of an image J[x,y), for / = 1, ID high-/ 
low-pass filtering is first applied to the original image 
f{x,y) along the horizontal direction (0°), followed by 
a decimation in which every odd-numbered element 
is removed. The ID filtering and decimation are then 
applied in the vertical direction (90°). According to the 
combination (high-high, high-low, low-high, and low- 
low filtering), the output consists of four quadrants (q). 
The lowest-resolution quadrant corresponds to the 
sub-bands for the scaling coefficient C±{k). The other 
three quadrants with directionalities {0°, 90°, 45"} are 
the sub-bands for the wavelet coefficients d\{k,q). The 
process is repeated for j = 2 by using Ci(/c) in place of 
the original image [13]. 

The coefficients d) (k, q) indicate the high-frequency 
elements of the decomposed image at position k and 
quadrant q within the frequency band /, where / is usu- 
ally referred to as the decomposition/ resolution level 
(or simply as the "level"). As the level / increases, the 
structural information of the image object in dj {k, q) de- 
creases [13]. Therefore, the wavelet coefficients used for 
the discrete WT may be chosen according to the size of 
the structure that one wants to highlight in the image. 
Alternatively, the images may be reconstructed to im- 
prove algorithm performance via MRA [14-16]. 



Algorithm optimization method 

Pixel intensities (in Hounsfield units, HU) of each slice 
were studied using MatLab 8 platform. The gray intensity 
levels of the pixels in regions containing enhanced and 
normal liver tissues were analyzed. The pixel intensity 
distribution in each type of tissue was fitted by Gauss- 
ians and the mean and SD determined in the slices was 
determined, as shown in Figure 1. The curve for normal 
liver tissue is depicted in part (a), contrasted liver tissue 
distributions in part (b), the superposition of (a) and 
(b) distribution is represented by the curve (c), and (d) is 
actual histogram extracted from the image. 

A virtual phantom was developed for the algorithm 
optimization. The phantom was used to optimize the 
detection of HCC and remarking the differences from 
normal liver tissue as shown in the phantom image of 
Figure 2(A). 

When constructing the phantom, the Gaussian distri- 
bution representing the normal tissue was used to fill a 
256 x 512 -pixel field image as a background. The pixels 
intensities for this environment were simulated according 
to the distribution represented by curve- (a) in Figure 1. A 
set of HCC lesions were simulated by cycles with max- 
imum diameter from 5 mm to 100 mm incrusted to that 
background. The circle areas were filled with pixels with 
intensity pseudorandomly generated by the Gaussian curve 
(b) in Figure 1. Several algorithms and wavelets filters for 
segmenting and quantifying the image structures were 
used until to get the best results. The final algorithm 
and filter configuration are described in next section. 

When calculating the efficiency of the algorithm, the 
diameters of the created circles in the liver were com- 
pared with the diameter measured by the algorithm. 
Circles of maximum diameters varying from 0.5 cm 
to 14.0 cm, in steps of 0.5 cm and 10 iterations each size 
were generated and used as input to the algorithm. This 
comparison is shown as a scatter plot located at Figure 3. 
Bland-Altman Limits of Agreement (LoA) encountered 
were in the range of -0.32 cm and 0.31 cm and an R 
squared equals to 0.99 for a linear fit, which is an accept- 
able limit of agreement. 

Optimized algorithm description 

The algorithm reads the DICOM image of the slice of 
interest, and the operator makes a manual segmentation of 
the liver tissue. For this step, a non automatic process is 
used, because the set of analyzed CT were previously seg- 
mented in the routine service. In fact, there are several al- 
gorithms available that can semi automatically segment the 
liver tissue with good results [17]. Of note, HCC is highly 
related to other liver diseases (i.e., hepatitis, venous throm- 
bosis, cirrhosis) that may deform the liver structure. Being 
so, we remark that manual segmentation of the liver tissue 
avoids the quantification of unwanted structures. 
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Hounsfield Units ( HU ) 

Figure 1 Distributions of pixel intensity for liver tissues. Distributions of pixel intensity for a real image containing normal curve-(a) and 
contrasted liver tissue curve-(b). the Gaussian distributions present mean value and dispersions 55 ±10 HU for normal tissue and 90 ±11 HU 
for contrasted liver tissue. The superposition of (a) and (b) distribution is represented by the curve (c), and (d) is actual histogram extracted 
from the image. 



The discrete WT frequency bands were estimated in 
ID by using equation (2) [18]: 



2 s p 



(2) 



where f s (^) is the center spatial frequency (mm 1 ) in 
the scale 5, p is the pixel pitch, and £o(f) is the 



pseudo-frequency (mm" 1 ) of the wavelet in its basic 
level (f 0l0 ). In this study, Daubechies 10 (£ 0 (^) = 0.693) 
[16,19] was used as the mother wavelet in all of the 
procedures. This asymmetric, orthogonal wavelet has 
many image-processing applications and shows better 
qualitative results than wavelets from other families 
[20-27]. The frequency bands chosen to reconstruct the 
ROI images were determined by considering the size of 
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Figure 2 Virtual phantom and algorithm results. Virtual phantom constituted by simulated liver tissue and encrusted carcinomas (circles with 
diameters of 10 cm, 8 cm, 6 cm, 4 cm, 2 cm and 0.5 cm) in part (A). This phantom was used to optimize the algorithm performance. In (B) an 
illustrative example of the algorithm performance without wavelet filtering and in (C) an illustration of the algorithm performance to highlight 
the HCC simulated area. 
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Figure 3 Algorithm results when applied to the phantom. Scatter plot of differences between algorithm measure of maximum diameter of 
the masses of the phantom (Figure 2), an agreement of ±02 cm is clearly shown. The central line corresponds to the mean value of deviations. 
The upper and lower lines depict the limits of 2 SDs. 



the HCC lesion. Approximation coefficients above level 
2 (= 0.22 mm) were used to reconstruct all of the im- 
ages, because such coefficients remove small-sized fluc- 
tuations (but not small-sized tumors) from the image. 

The output of the previous step was binarized with 
a threshold of M NO rmal + STD NO rmal» as stated in 
equation (3), 

Rl ( v v \ - I 1 if I ( X 'y) > M NORMAL ,~s 
B outputl*.W - jo else W 

where / (x, y) is the input image. After binarization, ero- 
sion (4) and dilation (5) filters were applied in the image 

B output(^^) : 

B 2 (xv) -!° ifBi utput ( X + i --2,y+;-2) = lforalli,j = 1,2,3 

(4) 

lW^M-j 0 else 

(5) 

In these equations, i and j are arbitrary labels for the 
neighbors of the pixel being analyzed. Erosion and dila- 
tion filters (when applied sequentially in an image) have 
the property of smoothing the objects in the binary 
image and removing small-sized objects produced by 
binarization [15,28]. 

Statistical assessment 

Three equivalent groups to the measured diameter by 
the radiologist were used: (i) The Gl, for which nodules 
equal or less than 3 cm in HCC diameter; (ii) The G2 
with nodules between 3 and 5 cm in HCC diameter 
and (iii) The G3, where HCC diameters were above 5 cm. 
A group of three experienced radiologists, with more 
than 5 years of experience each, scored the images to- 
gether agreeing with each other in the measure of the 
maximum diameter of the lesion. The separation 



criteria were based on the limits stipulated by Barce- 
lona Clinic Liver Cancer (BCLC) staging, in which 
multinodular lesions greater than 3.0 cm or unique le- 
sions greater than 5.0 cm are intermediate HCC and 
cannot be submitted to liver transplantation before 
tumor reduction [3]. Dependency was evaluated by the 
R-square produced by a linear fitting of computed 
against radiologist measures. T-student test was applied 
to conclude if there is statistical difference between the 
groups of measures and Bland-Altmann plots were 
also used to assess the dispersion and limits of agree- 
ment between measures [29]. 

Results and discussion 

The optimized algorithm was used to detect HCC in the 
liver of actual CT images. Figure 4 shows some examples 
of the algorithm performance for highlighting the HCC 
region. The algorithm located and dimensioned the le- 
sion in all images with high precision. The radiologist 
opinion was in quite good agreement with the delimita- 
tion region of the tumor for all the data set. 

Validation of the algorithm results was carried out by 
comparing measured diameters using the algorithm with 
the radiologist evaluations. Data were divided in three 
groups according to the maximum diameters calculated 
by the radiologist, as stated in the statistical assessment 
session. Figure 5 compares the results obtained from the 
algorithm computed measures with the results evaluated 
by the radiologist, showing the difference between the 
algorithm and the physician's opinions. The value of 
R-parameter around 0.97 was obtained for all groups. 
The expected linear dependency between the radiologist 
measures and the algorithm was observed. The T-student 
test did not find statistical differences between any group 
(p > 0.05). Figure 6 compares Bland-Altmann LoA for the 
three groups. It is clear that LoA tends to increase with the 
size of the tumor. 
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Table 1 shows mean and differences obtained by the 
radiologist and the algorithm. The algorithm measure- 
ments presented minimum errors in the group of 3 cm < 
d < 5 cm diameter. Larger HCCs presented greater necrotic 
parts inside the tumors, which makes its maximum diam- 
eter estimation more difficult by radiologist. 

Of interest, the algorithm was able to measure area 
and volume of the tumors evaluated. This can be an im- 
portant tool contributing to a better assessment of the 
tumor in the clinical routine. 




2 4 6 8 10 12 14 
Radiologist Measure of Maximum Diameter (cm) 

Figure 5 Radiologist and algorithm correlation in measures. 

Comparison between the maximum diameters determined by 
radiologist and by the algorithm. Diameter evaluated by radiologist 
versus diameter evaluated by the algorithm linear fit obtained an 
R-square of 0.97 and t-student test presented a p of 0.1 3, not indicating 
difference significant. 



Conclusions 

Due to the relevance of the tumor size to choose the 
best treatment for each patient with HCC lesions in the 
liver, the maximum diameter of the HCC lesions should 
be considered as a measurement that requires other re- 
sent complementary measurement methods. When the 
tumor maximum diameter is close to the limits estab- 
lished for the OLT inclusion criteria and there was not 
tumor spread outside the liver, TACE is one of the most 
treatment modality used to reduce HCC lesions, but the 
evaluation of contrast-enhanced exams after TACE proce- 
dures can be difficult because the lesions tend to be many 
necrotic areas within the tumor tissue. Since there are 
strict criteria used to indicate which treatment must be 
performed to patients with HCC and the most used criteria 
are based on the tumor size, a non-subjective evaluation is 
an important tool to be considered in this case. 

Currently, the new concept of tumor viable tissue is 
coming to improve the assessment of HCC lesions, espe- 
cially after TACE procedures. Some studies showed that 
mRECIST and EASL criteria are suitable to evaluate 
these tumors, being now considered more adequate than 
the previous WHO and RECIST. They are more closely 
correlated to clinical endpoints [30] and are easier to 
use. However, they can lead to variations in the evalua- 
tions obtained from different radiologists, whom are not 
able to calculate the viable tumor area. 

The analysis of the viable tumor area can be very import- 
ant after TACE, discriminating the lesions that really re- 
spond after the procedure. TACE often lead to coagulative 
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Figure 6 Algorithm bearing for small, medium and large HCC enhancement. Bland-Altmann plots of the nodules equal or less than 3 cm in 
HCC diameter (A), nodules between 3 and 5 cm in HCC diameter (B) diameters above 5 cm measured by the radiologist (C). The central 
lines corresponds to the mean value of deviations. The upper and lower lines depict the limits of 2 SDs. 



necrotic lesions in the central area of HCC tumors, and it 
can be difficult to be measured using only the current cri- 
teria [30]. For instance, a tumor of 5 cm with a necrotic 
are of 1 cm in the central area is different of a tumor 
of the same size with a central necrotic area of 4 cm, but 
both have the same maximum diameter according to the 
current criteria, that take into account the longest meas- 
ure of viable tumor tissue and not the necrotic area. On 
the other hand, the analysis of the viable tumor area can 
show the difference in these two tumors after a TACE 
procedure because it calculates the area of viable tumor 
in each lesion. Additionally, it can be used to compare 
the efficacy of TACE procedures using different chemo- 
therapeutic agents or of each hemodynamic service. 

Therefore, a single dimension does not take into ac- 
count the viable area of a given tumor, and better evalu- 
ation methods are needed to avoid misinterpretations 
about each nodule found at the contrast-enhanced exam. 
In this study, we have presented a novel method to meas- 
ure contrast-enhanced HCC nodules in liver CT exams. 
Our method was able to distinguish the normal liver tissue 
from the cancer, by using wavelets base to determine ac- 
curately the tumors' limits. This image treatment plays an 
important role in the lesions characterization, properly 
measuring the size of HCC nodules. The algorithm was 



able to measure small-sized tumors with great accuracy. 
The precision is gradually decreasing for large-sized le- 
sions (diameter > 5 cm), but still with enough precision to 
apply the used treatment criteria. 

Our results confirm the findings of previous studies, 
suggesting that the maximum diameter should not be used 
alone to represent the tumor size. Many HCC tumors are 
noncircular lesions and may be only approximated by an 
ellipse, as described by Jensen et al. [31]. Based on these 
findings, patients with different tumors' shapes can be mis- 
evaluated if only the maximum diameter would be used to 
classify the disease according to the mRECIST criteria. 
Therefore, the HCC tumors evaluation can be improved if 
the total viable tumor area could be analyzed, especially for 
tumors submitted to TACE procedures, which often have 
coagulative necrotic areas inside of them. This kind of 
evaluation of the local response after TACE could lead to a 
more detailed comparison between the local response and 
the clinical end points, but it needs to be evaluated 
in further studies. The trend to major difference for 
lesions >5 cm may be associated with necrosis caused by 
TACE, leading radiologists to overestimate lesions. For the 
moment, the algorithm presented in our study can be very 
useful to analyze the viable tumor area without subjective 
measurements and assist in the final decision clinic. 



Table 1 Difference between radiologist and algorithm estimation of maximum contrasted diameter of HCC 



Group 


Images 






M ± 2 x SD (cm) 


R a 




Radiologist 


Algorithn 


i Algorithm and radiologist difference 




d < 3 cm 


12 


2.2 ± 0.5 


2.1 ±0.6 


0.1 2 ±0.58 


0.78 


0.17 


3 cm < d < 5 cm 


28 


4.0 ± 0.6 


4.0 ± 0.7 


- 0.05 ± 0.86 


0.61 


0.54 


d > 6 cm 


23 


8.4 ± 2.7 


8.2 ± 2.4 


0.25 ± 1.14 


0.96 


0.05 



a R-Squared correlation coefficients between data and linear fit. 
b p values based on t-test. 

Difference between radiologist and algorithm estimation of maximum contrasted diameter of HCC. 
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